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We introduce a new family of integrable theories with N bosons and N freely ad- 
justable mass parameters. These theories restrict in particular limits to the "generalized 
supersymmetric" sine-Gordon models, as well as to the flavor anisotropic chiral Gross 
Neveu models (studied recently by N. Andrei and collaborators). The scattering theory 
involves scalar particles that are no bound states, and bears an intriguing resemblance 
wih the results of a sharp cut-off analysis of the Thirring model carried out by Korepin 
in (1980). Various physical applications are discussed. In particular, we demonstrate that 
our theories are the appropriate continuum limit of integrable quantum spin chains with 
mixtures of spins. 
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1. Introduction 

A variety of low dimensional experimental condensed matter systems have been stud- 
ied recently, that involve field theories with several bosons. Examples include tunneling in 
quantum wires, where two bosons are necessary to describe the charge and spin degrees 
of freedom of the electrons QXJ] , tunneling between multiple edges in fractional quantum 
Hall devices 0, nanotubes and two-leg ladders ||, etc. Properties of interest in these 
systems are usually non perturbative, and only a few techniques are available to obtain 
quantitatively reliable results, mostly conformal invariance and integrability. The search 
for integrable quantum field theories with several bosons is thus of some importance. 

The problem is, that besides the sine-Gordon model, most known integrable bosonic 
theories are of little practical use: they are usually of Toda type, and involve real exponen- 
tial of fields, that usually do not appear in a condensed matter context. Some exceptions 
to this unsatisfactory situation are known: for instance, the double sine-Gordon model 
turns out to be exactly solvable for some values of the couplings H , , , with potential 
applications to quantum wires. Also, the "generalized supersymmetric" extensions of the 
sine-Gordon model |7] can be rebosonized using standard bosonization formulas for the 
parafermions [§. These theories are useful in the discussion of the multichannel Kondo 
model |§, []T0|j , [|TT| ; the Af = 1 supersymmetric sine-Gordon model also appears in the 



context of quantum wires [|12] 



In this paper, we point out that there is a simple, integrable family of theories extend- 
ing the generalized supersymmetric sine-Gordon models, that involve iV bosons and have 
N adjustable mass parameters. This family can be considered as an extension of the flavor 
anisotropic Gross Neveu models that have been studied in the last few years by N.Andrei 
and collaborators ( mostly in the context of the channel anisotropic Kondo model |l3 ], [14] , 
||15|| ), to the case where an anisotropy is introduced both in the color and flavor sectors. 
The case where the color anisotropy is at the special "Toulouse" value is of special interest 



for applications to quantum wires or dissipative brownian motion |T6 



The models are presented in section 2, where integrability is proven and various limit- 
ing cases discussed. The scattering theory is discussed in section 3. The "classical" limit is 
analyzed in section 4, providing a general check of our approach. In section 5, the relation 
with quantum spin chains involving several species of spins is discussed. Some applications 
to impurity problems are discussed in section 6. Some final remarks are collected in sec- 
tion 7. The appendix contains numerous details on the numerical treatment both of the 
perturbation theory and of the TBA. 
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2. Generalities 

2.1. The integrable theories 

We consider a system of N chiral bosons with propagators 

N - 1 

< <Pj(zm{w) > = - 2 ^j— H z - w) 

2 

< (f)j(z)(f)k{w) > = — ln[z - w), j ^ k, 



(2.1) 



and introduce the following fields: 

= -L I ^V fe e^ fc ) , (2.2) 



xk = l 



'N 

where to = e 2l7T / N . These fields provide different realizations [§ of the fundamental 



parafermion |17| of Zjv type. The bosonic fields <fij are not independent (one can set 
indeed 4>n = ~4>i — ■ ■ ■ — 4>n-i)', they can be expressed in terms of N — 1 independent 
fields $j obeying 

< ^(z^H >= -25 jk ln(z - w), (2.3) 

by the transformation 

<j> j = e j9 $,j = l,...,N, (2.4) 

where the • denotes scalar product and the ej are weights of the fundamental representation 
of SU(N$. 

Introduce one additional bosonic field, which has trivial contractions with the preced- 
ing ones, and obeys 

< $(z)$(w) >= — — In(z-w). (2.5) 

4-7T 

Consider then the action (we assume all the aj are real positive numbers) 

AT-l 



A =i J dxdy £ [0„ + 2 + % (* + * 



\ / \ ( 2 - 6 ) 

+ ( ^a^(*) j I J2 a ^ e^ z)+ ^] + conjugate. 



2 That is, ei = Ai, ei = A2 — Ai, ejv-i = Aat_i — An-2, ejv = — Ajv-i, Ai the fundamental 
weights of St/ (iV) and <I> the (JV — 1) dimensional vector of coordinates $1, . . . , $jv— 1. 
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In the case where all the coefficients aj but one vanish, this is the action of the "generalized 
supersymmetric" sine-Gordon model, which is known to be integrable [[F]. We claim that 
this only a particular case of a more general integrable model, given by ( [2.6| ). 

To establish this result, we first observe that the fields obey the short distance 
expansions 



v[/0') 



where 



(z)¥ j \w) 



1 



2(N-1) 
[Z — W) N 



l + 2^-^(z-w) 2 T^\w) + . 



(2.7) 



1 



N 



N + 2 



k=i 



k^i 



and 



where 



w 



yz — w) n 

. N 



z-w)J (jk) (w) + . 



[2.1 



k=i 



We can then prove integrability, following , by establishing the existence of non local 
conserved currents. Introduce 



(2.9) 



then the short distance expansion of this current with the first term in the action reads, 
for the chiral part, 



E&j ] ( ^flfc* (j) (w)) exp (-i^--^$(z) ) exp (ip$(w)) 

J=l ) \k=l 



(3 N 



1 



TV 



(z — w) : 



^2b k a k + ( 



k=i 



z -w)^2 b k aiJ (kl) (w) + ...\ exp 

k^i 



8in 



$(z) +i/3$(w] 



The residue of the simple pole with thus be a total derivative iff the factor of (z — w) in 
the first bracket vanishes. This is equivalent to the condition 



£& fc o,[>-'> m - 1 
k^i 



0,to= 1,...,N- 1 
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which always has solutions, since it is a system of N — 1 equations with N unknown i. 

The short distance expansion of this current with the chiral part of the second term 
in the action has a leading term that goes as (z — w)~ 2 / N (z — w) 2 / N , and thus no simple 
pole. Following the standard argument, the expansion of J~ having a simple pole whose 
residue is a total derivative, the non local charge J J~ is conserved to first order in the 
perturbation. For generic value of (3, one can then argue that this is true to any order in 
perturbation theory, and, presumably, non perturbatively as well. 

Another conserved current is easily found by complex conjugation: 



J2b*^(z)\ exp (i ■ (2-10) 



The conservation of J then ensures integrability |18[ . 



2.2. The case N = 2 

Let us discuss in more details the simplest example where N = 2. In that case, the 
parafermions are self-conjugated (up to a sign). We set 

?/> < - 1 * ) = — i\[2 sin 0i = %x 

(2.11) 

-j/;( 2 ) = ^/2COS01 = 1p, 

where ijj and x are (real) Majorana fermions. The perturbative part of the action reads 
then 

{mj) + ib X ){a$ + ibxY m + {aif> - ibx)(a$ - ibx)e~ m \ (2-12) 
that is, regrouping terms 

2 [a 2 V>$ - b 2 xx] cos (3$ + 2ab (^X + x$) sin P®- ( 2 -13) 

The non local conserved currents read then 

J+(z) =(aijj - ib X ) exp (i^J®} 



4:77 

J (z) =(aip + ibx) exp ( — i— $ 



(2.14) 



3 The solution is easily expressed using the matrix N x (N — 1) matrix M whose elements are 
Mjk = aj-ffe_i by bj equal to the j th cofactor. 
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If a = or b = 0, the action reduces to the one of the supersymmetric sine-Gordon 
model (with an additional, decoupled, Majorana fermion). If a = b, the combinations 
appearing in the action become Dirac fermions, ip + %x = V^e 1 ^ 1 . We can thus reexponen- 
tiate them, to write the perturbing term as acos(/?$ + 4>i), so the model is equivalent to 

(0') 2 B 2 1 

a sine-Gordon model, at a coupling constant j3 with = + f ■ 
The currents have a fractional spin s = —, where 



In the case a = or b = 0, the currents are generators of the algebra sl(2) q |L8], with 
deformation parameter q = —e~ l7r /' r . In the general case however, they do not form a 
closed algebra. 

To understand the situation a little better, it is useful to go to the SU (2) symmetric 
point f3 2 = 4tc. There are two underlying level one algebras, with generators 



4> + iX 
V2 



-e' 



= i^rx e - t ^ (216) 
y/2 



jj 5 =i (-tpx + V^Ttd^j 



and 

J- =t±3 e -^ (2.17) 
y/2 

and all short distance expansions between operators of different algebras are non singular. 
The sum J\ + J<i provides a level two representation. For general a, 6, the currents can be 
written as combinations of the J\ and J2. Setting a = fi + A and b = fx — A, we have 

J+ =\J+ + 

J~ =AJ 1 ~ + 11J2 ■ 

This is suggestive of a system where two flavors of fermionic currents are combined in a 
flavor anisotropic fashion. Indeed, introduce new bosons defined by 

(2.18) 

(fil + V47T$ =(f>2, 
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the perturbing term is then proportional to 



A 2 cos(v?i + ipi) + fj 2 cos(v?2 + V2) + A// cos((/?i + ^2) + A/i cos(y?i + ip 2 ) 

This is the abelian bosonized form of a Gross Neveu type interaction with flavor anisotropy 

(XJf + //Jf ) x (XJf + //Jf ) + (x^y) 

The zz term is missing in this interaction - it is well known that this term is generated 
under renormalization |J. Away from the SU(2) point, one can similarly consider our 
model as a color and flavor anisotropic chiral Gross Neveu model (upon bosonization, this 
model gives rise to 4 independent fields, but only 2 of them appear in the interaction due 
to chirality, the other 2 ones contributing free parts to the action). 

2. 3. Arbitrary N 

The previous discussison easily extends to other values of N. The currents have 
a fractional spin s = -, where 

7=^V- (2-19) 

N 8tt 

In the case a = or b = 0, the currents are generators of the algebra sl(2) q with deformation 
parameter q = — e -Z7F / 7 . In the general case, they do not form a closed algebra. In the 
limit [3 2 = 47r, they can be expressed as combinations of iV generators belonging to iV 
different realizations of a level 1 SU(2) algebra. The rebosonized action is a Gross Neveu 
models with two colors and N flavors, and flavor anisotropy, with an interaction term of 
the form 

where the anisotropy coefficients Xj are related with the terms in the original action by 



Xj = ^u 3k a k . (2.20) 



k=i 

That the Gross Neveu model with flavor anisotropy is integrable has been pointed out 



several years ago in fact, in [^3[,[^3[. Integrability is established there by direct diagonal- 
ization of the bare hamiltonian together with "dynamical fusion" . An intriguing feature is 
that the proof presented in WM , strictly speaking, works only for the case N = 2 (and some 
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subcases of special flavor anisotropy for larger N). The reason is, that in the approach of 
I|J, the bare particles must have a bare flavor scattering matrix that is a solution of the 



Yang Baxter equation. In the flavor isotropic case, this S matrix is the standard SU(N) 
R-matrix; but, as far as we know, there is no way to deform this R matrix in a non trivial 
way by introducing N — 1 independent anisotropic parameters - all available single and 
multiparameters quantum group approaches still explore a very small subset of all the pos- 
sible flavor anisotropies. On the other hand, from the point of view we have adopted (that 
deals directly with the renormalized action), all flavor anisotropies play equivalent roles, 
and integrability appears generally true. The argument also extends straightforwardly to 



the case of color anisotropy, not considered in |T3 



Of course, a particular choice of anisotropy is when one of the coefficients Xj vanishes 
exactly, in which case the modele reduces to one with N — 1 flavors. In the case of general 
/3, the same conditon, say A/v = leads, using formula (pT2|), to a problem where the field 
(f>N has disappeared from the action. We are then left with a set of N — 1 fields satisfying 

N - 1 

< <i>j{z)(t>j{w) > = -2 — ln(z - w) 

N -2 2 

2 1 ' (2-21) 

< (j>j(z)(j> k (w) > = — ln(z - w) 

° 2 

\n{z — w) — -Trr-Tr rr hi(z — w),j ^ k 



N-l v ' N(N-l) 

We can then write <pj = <// + where there are N — 1 fields satisfying relations similar 
to ( |2.21| ) but with the replacement N — > N — 1. The problem is then equivalent to the 
case N —y N — 1, but with a shift of the leftover exponential, § > § — I 



3. Conjectured scattering theory 

3.1. The thermodynamic Bethe ansatz 

Though there are arguments based on symmetry to infer what the scattering theory 
should look like, the approach we use is to first conjecture a set of thermodynamic Bethe 
ansatz equations (TBA) to compute the free energy of the 1+1 quantum field theory 
associated with the action ( |2.6| ) at temperature T. We parametrize 



and consider the case 7 an integer. Our conjecture is as follows. Introduce the TBA 
diagram 

O -y + N-1 
1 2 N N + l J 

<S> <8> O O O 7 + N-2 

O l + N 



with incidence matrix Nj k such that Njk = 1 if the nodes j and k are connected, and 
otherwise (in particular Njj = 0). With this diagram, we associate the set of pseudo 
energies (one for each node) solution of the system (R — 1/T, T the temperature) 

«,=£) * jWn »Hc»h»-5>, / ^-^1 ta(l + e -«C)). (3.2) 

fe=l fc V ^ 

The free energy reads then 

F = -— V / ^-m fc cosh01n(l + e- efcW ) . (3.3) 
27T * — ' I 2tv \ J 

fe=l J 

In the foregoing equations, the are a set of masses which depend on the couplings 
dk in the bare action. By dimensional analysis, [ak] = [length] 4^ ~ w . Therefore, we have 

JV + -Y 

m k = G k (a 1 ,...,a N ) = a x 2 F k (a 2 /a 1 , . . . , a N /a{). (3.4) 

The Gk are homogeneous functions of the couplings a k - Some properties of these functions 
are known before hand of course. They are symmetric functions of their arguments. If 
all the afc but one vanish, we know that the problem becomes equivalent to the N th 
supersymmetric sine-Gordon model, and therefore, from known results |[L8|1 , [[19]], all the 
masses but the N th one must vanish. This means that Gj = 0, j = 1, . . . , N — 1, when 
all the a k but one vanish. Also, we know that the N th mass vanishes when one of the 
fields decouples, ie when one of the coefficients (|2.20|) vanishes. More generally, the masses 
m k . . . mjv vanish when N — k + 1 of these coefficients vanish. We will get back to the 
determination of the functions G k below. 

The evidence for the TBA comes first from the compatibility with all limiting cases. 
Moreover, the analysis of the Y system [|19| associated with it shows that the dimension 
of the UV perturbing operator is always the same as in the generalized supersymmetric 
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case (it does not depend on the number of massive nodes), ie h = 1 t f ^ 1 = §^ + 
as desired. Finally, the central charge, in the generic case when all the mk are non zero 
is simply equal to the number of massive nodes, ie c = N. This is easily checked. Using 
standard formulas, the central charge is expressed in terms of the solutions of the system 
( |Q|) as T — ► and T — ► oo. In the first case, the iV first e's are all infinite, the others 
follow from 

Xj = = (j + l) 2 - 1, j = N + 1, . . . , TV + 7 - 2 
XN+-T-1 = xn+j =7 — 1. 
In the second case, one has to solve the same system with more nodes, ie 

Vi = e " £j = 0' + 1) 2 - i,i = i, • • • , ^ + 7 - 2 

y N + 1 -i = = N + 7 - 1. 
The central charge is then (here L designates the Euler dilogarithm) Jl9j 

For a D diagram, as T — > cxd, the sum ^ L f i+y) ^ s e Q ua l to the number of nodes 

2 

minus one, so the central charge is simply, from L(l) = equal to the number of massive 
nodes, ie c = N indeed. 

3.2. Scattering theory 

The scattering theory associated with this TBA is very simple. One introduces a 
set of iV — 1 scalar massive particles with masses mi, . . . , mjv-i- One also introduces a 
pair soliton/antisoliton with masses uin- The latter scatter with the usual sine-Gordon S 
matrix that corresponds to the quantum group parameter introduced above, q = — e _l7r / 7 
- it is the same as the S matrix of an ordinary sine-Gordon model at coupling — 



8tt 7+1 - 

The scalar particle of label k scatters trivially with all particles, except the ones of label 
k ± 1, with which it scatters with the CDD factor S = itanh (| — When k = N — 1, 
the particle scatters with the soliton and antisolitons with the same CDD factor. It is 
important to stress that the sine-Gordon S-matrix considered here has no poles in the 
physical strip: the scalar particles are not bound states of the soliton and antisoliton. 

Remarkably, a scattering theory built with similar ingredients appears in a paper by 
Korepin |2(J . There, the author discusses the Thirring model in the repulsive regime, using 
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a sharp cut-off regularization. For a coupling corresponding to a sine-Gordon parameter 
^ E -A-j-, l fp2 ■> ne finds, in addition to the soliton and antisoliton, a spectrum made 
of (/ — 1) neutral particles, with same S-matrices as ours, but where all the masses are 
uniquely determined as a function of I (in particular, the soliton mass becomes infinite 
when ^ = ). The soliton and antisoliton scatter through a sine-gordon S matrix with 
a renormalized J = l _ L ^ l ^i^f_^yl n ^ parameter. 

The relation with our problem, if any, is not clear. It is usually admitted that the in 
the repulsive regime of SG, the quantum theory must be defined with care, and depends 
on the cut-off. For smoother cut-offs, as well as XXZ type regularizations, the results of 



P0(1 are not supposed to hold, and the standard description of fl2]J with only the soliton 
and antisoliton to be correct. 

Observe that the sine- Gordon part of the S matrix commutes with the quantum alge- 
bra sl(2) . While the non local conserved charges Q ± have commutation relations that do 



not close, they can be expresed as combinations of iV basic charges generating sl(2) . A 
simple representation of the algebra generated by the is then obtained by identifying 
all these charges - it is realized on the multiparticle soliton antisoliton states as in the 



usual sine-Gordon model ||18| . In the general case, there are no other conserved charges. 
If N = 2 for instance, another, local, conserved charge appears when eg b = 0, since the 
local currents Q = -i/><9$ is conserved j!8| . But away from this value (and a — of course) , 



this is not true. We must thus complete the S matrix by a sector with no apparent sym- 
metries, that does not spoil the Yang Baxter equation: besides "vertex" and "RSOS" type 
solutions, the only available choice is a set of particles with diagonal scattering. Requiring 
the central charge to be N, and the TBA to restrict to the known ones when some of 
the couplings vanish (and additional symmetries appear) seems to leave no choice but our 
result. 

Still another check comes from discussing the limit /3 — > 0, to which we turn now. 



4. The "classical" limit. 

We call here classical the limit where the sine-Gordon part of the scattering theory 
becomes identical with the classical one, that is (3 — > 0, 7 — ► 0. 
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4- 1- The TBA at the reflectionless points 

So far we wrote a TBA in the simplest case 7 an integer, for which |^ > N ^ N+1 ^ , 
corresponding to the SG component of the S-matrix being in the repulsive regime. To 
approach (3 = 0, we consider instead the cases 7 = ^, that is ^ = jv(reAr+i) : the SG 

5- matrix is now in the attractive regime, at the so called reflectionless points. This means 
that the spectrum has to be completed by the bound states of solitons and antisolitons, 
the n — 1 breathers of masses 2m sin (jj^), where m is the soliton mass. We denote 
2r»i sin (^-) , . . . , 2rriN-i sin (^-) the masses of the scalar particles (recall these are not 
bound states; the factor 2 sin (^-) in their masses is introduced for convenience only). By 
building the complete scattering theory using bootstrap and fusion, one finds the following 
TBA equations. 

One first has equations for the right hand side of the diagram, that look like the 
standard ones for the attractive regime of sine- Gordon 

€j = N jk K * In (l + e £k ) ,N +1 < j < N + n. (4.1) 

k 

There is then a central part involving the nodes — 1 and A, with 

e N = K* [ln(l + e £jv+1 ) - In (l + e"^" 1 )] , (4.2) 



and 

/ tv \ ] / TV \ cosh 9 
2mN-i cos y— j + m tan ( — 



e/v-i — 



2n/ T (4.3) 
A * In (1 + e ejv ) — K" * (l + e"^" 1 ) — K' * (l + e~ eN - 2 ). 



Finally, the left hand side of the diagram looks in turn like the usual repulsive one 

e, = 2m, sin {£) ^ * (l + e"*) ,j<N — 2. (4.4) 

k 

These equations can be conveniently encoded in the diagram 



N + n-1 

12 A — 1 A 

0-/-0- - ®-h®- -0- ® ® A + n-2 

N + n 
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The asymptotic conditions for the attractive part are 



6j ~ 2m sin 



U-N + 1)tt 



2n 
cosh 9 



cosh# 



N < j <n+ N -2 



(4.5) 



en+N-i ~ e n+N re- 



introducing the fourier transform g(6>) = f °° l^e ^^(u;), one has -KT = 75 — — t^-v 
(the standard attractive kernel), = 2 cosh(>^/2) ^ ne standard repulsive kernel), K" = 

cosh(n— l)nui/2r 



2 cosh nui /2n cosh 7rtj / 2 " 

For T — > 0, we have Xj = 0, since now all nodes are massive. For T — > oo, we have 

=(j + l) 2 -l,j<N-l 

2 



j" - N + 1 + — ] - 1. N < j < N + /? - 2 



e £N+ n -l =e tN+ n =n ^ I 



1 

TV 



The central charge is thus 



(4.6) 



C = —r < 2L 



7Y- 



n + 



N 



n-1 



5=1 



& + J 



JV-1 
3=1 



(J + 1) 2 J 



iV. 



(4.7) 



^.2. The classical limit 

The "classical limit" is obtained by letting (3 — > 0. In our TBA, this means n — > oo. 
To get non trivial results then, we scale the soliton mass m with n, so the mass of the first 
breather remains finite. Similarly, we assume that the parameters rrij are also scaled with 
7i, so that the masses of the scalar particles remain finite. Our computation follows the 
general strategy of |]22| , | ]23|| . In that limit, it is convenient to introduce the new notation 
Kj = ej+N-i, j > 1. When n — ► oo, the kernels K and K" become delta functions, (we 
set K' = s) and one finds the general solution 



e Kj + 1 



aA j - a~ 1 A~ j 
A- A- 1 



The constant A follows from the knowledge of mass terms, A = e c ^ 2T ,C = ^cosh^, 
while the constant a depends on e;v_i: 



1 + e 



-ejv-i 



A — A' 



a — a 
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In that limit, the equation satisfied by e^-i is 



e^.! = -s * In (1 + e"^ 2 ) - 1 In (1 + e Kl ) - 1 In (l + e"^- 1 ) + (a + 



where A = 

Let us then define 



s*ln(l + e^-* ) + A^. 



By simple algebra, one finds that the following holds 

1 + e - ^- 1 = (1 + e" ejv ) (l + e" e ' w -i) 

together with 

C 

e N = (A + 1)- - s * In (1 + e" £jv - 2 ) 

We can thus trade completely the right hand side of the diagram for an additional 
node, and get the equations (recall s = K' ,s = 2c osh(7r^/2) ) 



_ M, cosh 9 
* = T 



where the TBA diagram is a D diagram with A?" nodes 

(g) N-l 
1 2 AT-3 / 

-0- -0- N — 2 

\ N 

and masses, in that limit are Mi = wii^, . . . , Mjv_i = mAr_i^, and M/v = Afjv_i + m^. 
We can now compute the free energy. Its general expression is 

? = ~t if 'i c m '« (i + «""') ~ £ / "» f 1 + ■ «•») 

By using the basic TBA equation 

Ki « 4r + 2^/cln(l + e" Kfe ) - In (l + e"^- 1 ) , 



fc=i 



13 



one finds after a few simple manipulations 

^ = -53/ ^M iC osh01n(l + e- e '-) + J cosh(^) In (l - e" M cosh9 / T ) . (4.10) 

Here, the new mass = Mjv_i + ^ as before, and M = = . 

The results of the classical limit are therefore equations ( |4.8| ) and Q4.10Q . This means, 
the classical limit is made up of N particles scattering in a non trivial way, plus a decoupled 
free boson. 

In the particular case of the generalized supersymmetric sine-Gordon model, all the 
masses Mj but the N th one vanish: the system ( |4.8| ) reproduces the well known TBA for 
Z/v field theories perturbed by the parafermion field |23}| , as is expected from letting (3 — > 
in the action. 

In the particular case N = 2, there is no node N — 2, and the TBA system is trivial. 
The system decouples into two free fermions of masses and (mi + m) — , and a free 
boson of mass This is again expected from the action, and a rather non trivial check 
from the point of view of the TBA. The mass for the free boson certainly arises from 
counter terms analogous to ones arising in the M = 1 supersymmmetric action (for 
more discussion about this, see next section). It also follows that, in the classical limit, the 
correspondence between masses in the TBA and bare couplings goes, assuming a 2 < o 2 , as 

mi oca 2 

(4.11) 

m oc(6 2 — a 2 ), 

where m is the mass of the first breather, and mi the mass of the scalar particle. In 
the more general case, the decoupled free boson presumably gets its mass from a counter 
term still analogous to what happens in the supersymmetric case ||25|| . The rest of the 
TBA corresponds to a non trivial theory, with N species of Zjq parafermions interacting. 
It is interesting to discuss in more details the case N = 3. There, the TBA is based on the 
diagram 

12 3 

-0- -0 

with central charge c = 2. Observe that this value can be obtained by 2 = | + + | , 
corresponding to the sum of the central charges for the Ising, tricritical and tetracritical 
Ising model (or, alternatively, the Ising and tricritical Ising models, and the 3 state Potts 
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model). Remarkably, the weight h = | of our three parafermions \E r< -- 7 - ) can be recovered by 
using fields of these minimal conformal field theories. The three fields 5 , ^iJ 10 ^^J^ 
and $S3 /2) $33 /10) $g /5) (where lower labels are labels of the Kac table, upper labels are 
the central chagres) do have conformal weights h = | . One might be tempted to infer that 
the foregoing TBA also describes the perturbation of this product of three theories by a 
combination of these three fields. This cannot be true however: it is easy to check that 
the operator algebra of these three fields cannot be reproduced by using our parafermion 
fields only, due for instance to the appearance of powers 1/15 and 19/15. More generally, 
a TBA like 

1 2 M 

-0- 

6 



has a central charge that can be written as c = c\ +C2 + ■ ■ ■ cm , where cm = 1 — (m+2)(m+3) 

VI- 
M 



The conformal weight of the perturbing operator is h = M M X , which can be reproduced 



by M fields of the form ${ c 3 m) $ 3 c 3 m+l) . . . , m = 1, . . . , M. While it is tempting to 
speculate that ther TBA does describe the product of M minimal models perturbed by a 
combination of these fields, this result does not seem to be true. 



5. Numerical check 

To check the validity of the TBA besides the qualitative features we just discussed, 
one needs to compare the result for the free energy (|3.3| ) to perturbative computations. 
We restrict here to the simplest case N = 2. Because the action (|2.13| ) has two free 



parameters, reflected in the existence of the two masses m\ and m.2, a full consistence 
check requires at least going to the 6 th order (odd orders vanish) in perturbation - a really 
complicated task, as discussed in more details in the appendix. The second order does not 
give any check but fixes a global scale. The fourth order does contain some information: 
consistency determines uniquely the relation between the masses and a, b. In fact, it is 
not obvious a priori that this relation will be physical: finding it involves solving some 
quadratic equations whose solutions might well be complex, establishing, in fact, that the 
TBA is not the right one. We have however always found solutions that are physical, 
indicating at least that the TBA is consistent to that order. Moreover, the general shape 
of the functions G we obtain can be argued to be the right one based on limiting cases, 
giving us some confidence in the TBA indeed. 
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As an example let us consider the first non trivial attractive reflectionless case, which 
corresponds to |- = ^ . We have numerically solved the TB A equations ( flip , ( [4.2| ) , ( fOl ) 
with iV = 2, n = 2 for different values of the mass ratio m/mi. The coefficients of the 
analytic expansion of the running central charge 



C(r, m/mi) = 2 + Br 2 + o^r 



k(2-g) 



(5.1) 



k=i 



are determined by fitting. As a check, the bulk terms of the (extrapolated) limiting cases 
m/mi = and mi/m = 0, i.e. the sine-Gordon and the supersymmetric sine-Gordon 
points, agree with the exact values B = 3/vT and B = within an accuracy of 0.1%. In 
Fig.l we give the result for the adimensional ratio of coefficients 

2 



12 



r(i/2 + 



8ti- ' 



r(i/2- 



r(- 



§1- 

4-7T ' 



r(i + 



4tt . 



0,4 



(5.2) 



function of the mass ratio m/m 



l- 



20 



15 



10 



10 



15 



20 25 



m/ml 



Fig. 1: TBA result for the ratio I 

This curve has to be compared with the same universal ratio determined in the ap- 
pendix with perturbation theory as a function of x = (k-/k + ) 2 , shown in Fig. 2. 

The limiting cases are the sine-Gordon model at x = and the supersymmetric sine- 
Gordon model at x = 1. 
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Fig. 2: Perturbative result for ratio I 



1.2 
1 

0.8 
0.6 
0.4 
0.2 



0.1 



0.2 0.3 0.4 0.5 



. 6 



m/ml 



Fig. 3: "Quasi- classical" behaviour of the ratio of coupling contants 

The value of the universal ratio / at the minimum point found with perturbation 
theory is = -2.374 in good agreement with the TBA value I^f = -2.366. Another 
important check is the value of the universal ratio at the supersymmetric sine-Gordon 
point: the perturbative result is Fgg y = 29.45 while the TBA value is lJ u B s £ = 29.49, 
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confirming therefore the standard analysis!!! of this point 

By solving the second order equation in (k-/k+) 2 one can extract the dependence 
of the coupling constants on the mass ratio. In Fig. 3 we show that the quantity R = 
(b 2 — a 2 )/a 2 is still, for our value P 2 /8tt = 1/10, almost a linear function of m/mi, like 
in the classical limit (3 2 /8n — ► where R = moreover, the slope is very close to its 
classical value of 1/2. 



6. Relations with multispin integrable lattice models 



The conjectured TBA appears very naturally in an a priori different context: the 
study of inhomogeneous integrable lattice models of XXZ type with a mixture of different 



representations. To explain this in a concise manner, we refer the reader to [26|, and use 
similar notations (though the matter is quite standard). We consider thus an integrable 
model based on sl qo (2) R-matrices, whose "vertical space" is an array with spins s± = 
S2 = 1, S3 = S4 = 2, . . . , S2N-1 = $2N = N and S2N+1 = s« otherwise, ie made of 
blocks representing the iV first values of 577(2) spin. The associated spectral parameters 
alternate u\ = —U2 = i A ~t, Xl , ■ ■ ■ , M2JV-1 = ~U2N = i A+ 9 Ajv , and U2N+i = Ui otherwise. 



The anisotropy is determined by the quantum group parameter q Q 
hamiltonian we chose 



exp 



7+^ 



For 



H 



-ld_ 
t du 



\Tvlt' 



A + A x 
1 h u 



v . A + A;v 
. . . t ( 1 h w 



.iV 



,A + Ai 



.A + Xn 
-i u 



u 



(6.1) 



u=0 



Here, t s denotes the transfer matrix based on the foregoing vertical space and a "horizontal 
space" is a representation of spin s. The whole geometry can be illustrated on the following 
picture 



ui . . . 



u 



3 • ■ ■ 



U 



Si... 



J J . . . 



That is, no counter term is necessary to make the action supersymmetric and integrable away 



from (3 = 0. 
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In the case of an array with a single type of spin j, the physical equations are well 

ft 2 

known to be encoded in a TBA identical to what we studied above for |^ = N ^ +N ^ , but 
with a mass term on node j only. In this more general case, it is straightforward to show 
that one gets now a mass term for each of the N first nodes, with masses 

m^Mexp^-M^, (6.2) 

where M = | exp a) . 

By taking the continuum limit, this local integrable lattice model will give rise to an 
integrable quantum field theory that has exactly the TBA conjectured in section 3. This 
indicates that this TBA is more than an abstract set of equations, but must be related to 
a genuine quantum field theory: we conjecture this field theory is nothing but ( |2.6|) . 

It is interesting to observe that the same TBA would also be obtained by chosing 
a uniform spectral parameter (eg all \j = 0) but by putting different amounts of the 
various spins, with densities proportional to the masses rrij] see |27| for more details on 
this approach. In the latter reference in particular, S matrices are directly derived from 
the lattice regularization, and agree with the results of section 3. In the context of lattice 



models, central charges have also been computed with TBA similar to ours [23 , p9| 



7. Impurity problems 

The same argument of perturbative integrability carries through in the case of impurity 
problems | 3~0|] . Consider thus the problem with free bosons in the bulk, and an 



interaction term at the boundary 



N 

H bdr = ( ^a^ W (0) ] S~e l ^ 0) + conjugate. (7.1) 

Here, are raising and lowering operators in a spin representation of the quantum 
group sl(2) qo (here the deformation parameter is not the same than the one appearing in 
the S-matrices earlier, but rather the one of the lattice model above, qo = exp ^ ^ N j ) ■ 
Note that only the right moving part of the fields appears in the action, but that the bosons 
0, $ all have Neumann boundary conditions, ie their right and left moving components are 



In our conventions, the fundamental has spin one. 
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identical at the boundary. Alternatively, one could thus express the boundary perturbation 
with the total fields, and exponentials of half the argument. A particular case is where 
the boundary spin is in a "cyclic" representation, and can be gauged away pi| , p2| . One 
finds then 

N 

Hbdr = expi {(f) j + + conjugate, (7-2) 

i=i 

with Xj defined in ( 2.20| ). 

In the simplest case of the repulsive regime, and for 7 an integer, the boundary free 
energy for spn j with j < 7 + N — 2 reads, as in the usual anisotropic Kondo problem , 

d9 1 



bdr 



2tt cosh(# - B ) 



ln(l + e _e '). (7.3) 



Here, the ej are obtained by solving the TBA system ( |3.2| ) in the massless limit; this 
means, one sends the masses to zero and the rapidities to 00, such that only right moving 
particles with dispersion relation e = p remain. In the TBA, the source terms are obtained 
by the simple substitution cosh# — > \e 6 . The "masses" are simply parameters with the 
physical dimension of an inverse length. The rapidity 6 b is such that m\e B oc a 1 2 . For 
the case of the cyclic spin, the free energy reads as (Ol) but with j = 7 + N. 

It is e S pecia„ y inte reSt in g to consider the free e" y in the oLioaJ Ifrnit 1 . There, 
the same formula (|7.3| ) holds for spin j < N — 2 and the TBA now given by ( f4.8| ). For the 
spin j = N — 1, we have, due to some of the foregoing changes of variables, 



bdr 



Finally, for the case of cyclic boundary spin. 

F bdr = -T [ 1 ln(l + e-^), (7.5) 

J 2tt cosh(u — v B ) 

(in the last formulas, we have subtracted the trivial boundary free energy of the free boson 



In the case N = 2, the problem has been studied in ||34j| . Take an impurity of spin 
1/2, and set S + = ci + , S~ = d. Introducing the Dirac fermion T = i/j + the boundary 
action reads 

^— - (T+d+ + Td) + ^ (Td+ + T+d) 



This is usually called the "Toulouse" limit in the context of impurity problems |3 
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According to (|7.4[) , the free energy in the spin 1/2 case (ie j = 1 in our notations), reads, 
using notations of section 2, 



bdr 



d6_ 
2^ 



1 



+ 



1 



_cosh(#-# B ) cosh(0-0^ 



BJ 



In 1 + e" 



'/T 



where we have made a shift of the variable of integration, and thus e® B /e® 3 oc In 
terms of the A, \x variables describing the couplings to the different channels, this reads 
e °B j e 9 B f^jjl , in agreement with results of fl34fl . 

For general N, and when all the masses mk but myv vanish ( the standard gener- 



alized supersymmetric case) we obtain from ( [7.5| ) the ratio of degeneracy factors in the 



UV and IR 



guv 
gin 



N, corresponding to a flow from free to fixed boundary conditions 



in the Zn model. In general, we have here the solution of a multiboson problem with 
arbitrary couplings; when N = 3 for instance, this means that the problem with boundary 
perturbation 

A cos 4>i + fx (cos fa + COS (f>s) 

is integrable. Rexpressing the bosons (p m terms of the independent bosons $, we obtain a 
two boson problem that is nothing but a quantum wire problem: see |TB[ for more details. 



8. Conclusions 

We feel there is more to understand in the theories we have addressed. The relation 



with the bare Bethe ansatz solutions of []T3|], JTJ] in the color isotropic case is poorly 
understood, as discussed in the text. The scattering theory we have proposed is rather 
mysterious, and we have not answered the question of what the scalar particles have 
to do with the flavor symmetry breaking. The problem of analytically determining the 
relation between the action parameters aj and the masses rnj remains in general open. 
Finally, we have restricted to positive coefficients a,j in the problem, but it is clear that 
the perturbations will not always be massive if we allow some of these coefficients to be 
negative: where the theories flow to in that case is also an open question. 

Acknowledgments: We thank N. Andrei F. Lesage, and Al. Zamolodchikov for many 
useful discussions. 



21 



9. Appendix 

As a check of the correctness of our solution we now compute with conformal pertur- 
bation theory the free energy for the case N = 2. The result has to be compared with 
the free energy obtained by numerically solving the TBA equations. Let's consider the 
system in the strip geometry (R,L) defined by the action A = A c f t — f strip <& int , where 
the interaction is given by (|2.6| ). The dimensionless running central charge 

6 (it, a, b) 
becomes in perturbation theory 



lim — In Z[R,L] 



n 



(9.1) 



C(R,a,b) 



R 



c uv + 12 lim 

k=2 



U 



d 2 w\ 



■ d 2 W k < <S> int {W X ) ■ ■ ■ $ in t(Wk) >' 



strip 



(9.2) 

The first correction k = 2 is ultraviolet divergent for any real value of /3 2 , the anomalous 
dimension of the perturbing operator being A = 1/2 + Whatever regularization we 



choose, for example a radially ordered one P5| , we obtain a divergent part, to be subtracted 
by a constant counterterm in the lagrangian, and a universal part. The counterterm 
contains a possible finite contribution giving rise to a non-universal bulk term which has 
to be fixed by a normalization condition. This condition is C(R, a, b) — > when R — > oo. 
In practice the bulk term will be determined by comparison with the TBA result, while in 
integrable theories with only one mass scale it can be computed analytically. Taking into 
account the first non trivial correction, the running central charge reads 



C{R, a, b) =cuv + Cbuik (R, a, b) + 

9 2 n 2 



6(2tt)^ 



r(l/2 + 



01 

8ti- - 



r(V2-g 



rf- 



4tt - 



r(i + E 



(a 2 + b 2 )R x - 



4-7T 



1 2 



(9.3) 



The theory depends on the pair of massive coupling constants (a, 6), or alternatively on 
the two masses (ml,m). Since we are going to compare our perturbative expansion with 
the TBA result it is useful to introduce the dimensionless coupling constants (/c + ,/c_) 
and to make explicit the dependence on the mass ratio: a 2 + b 2 = k + (m/rrii) m x ~^ and 



a 2 -b 2 



k-(m/mi) m 



2 

± 4tt 



Defining the dimensionless quantity r = mR, the second 



order running central charge becomes 

C(r, m/rai) =2 + c huik (m/ mi ) r 2 

r(i/2 + 



+ 



6(2ttY 



PL- 

8ti- - 



-i 2 



r(i/2 



8tt- 



471" - 



r(i + 



A-k> 



[/c_l_(m/mi)] 2 r 2 27r 



(9.4) 
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Let's consider the case f3 2 < 2n, which includes the attractive regime. The only 
UV divergence of perturbation theory is the one that occurs at second order. All the 
other perturbative contributions are UV finite. The third order correction is zero because 
the unperturbed three points correlation function of the interaction is zero by charge 
neutrality. In order to compute the fourth order correction we map the w-strip onto the 
z-plane, z = exp(i2irw / R) , and we express by Wick theorem the unperturbed four points 
correlation function 



G 



(a 2 -b 2 )' 2 (a 2 + b 2 ) 2 (a 2 + b 2 ) 2 



+ 



Z12Z34 Z 13 Z 2 4 Z14Z23 

to which we have to subtract the disconnected term 





Z\l\ 


234 




\Z13 




214 1 


^23 


224 



Si 

2tt 



+ Perm. (9.5) 



2\4 



D = 4(a 2 + b 2 ) 



212 234 



+ 



213 224 



+ 



214 223 



(9.6) 



It is useful to rewrite the correlation function as the sum of two pieces 



G 1 
and 

Go 



(a 2 - b 2 ) 4 (a 2 + b 2 ) 4 (a 2 + 6 2 ) 4 



Pl2 P34 



+ 



4\2 



F13 P24 



+ 



PU T 2 2 3 





212I 


|234 




I213I 


214 


223 




224 



Si 

2tt 



+ Perm. (9.7) 



(a 4 - b 
Z12Z34Z13Z24 
+ Perm. 



+ 



(a 4 



(,4\2 



2\4 



(a 2 + b 2 ) 



2l22342 14 2:23 2i3Z242 14 ^ 2 3 



+ c.c. 





212 




234 




\Z13 




214 




223 




224 



Si 

2tt 



(9.8) 

and to compute the perturbation integral as the sum of the two corresponding integrals 
f (G\ — D) and f G 2 , both separately UV- finite. As a result, the fourth order term is the 
sum C< 4 ) = C{ 4) + Cf } where 



C{ 4) = 12 lim -5- / d 2 Wl ■ ■ ■ d 2 w± [G 1 - D] 

L^oo 2lXL J strip 



strip 



and 



= 12 lim — 

L^oo 2nL 



d 2 wi ■ ■ ■ d 2 W4 G2 



strip 



strip ' 



(9.9) 



(9.10) 



In the computation of the first integral we take the limit L — > 00 first, and we cancel 
the overall volume of the strip RL, and then we map the infinite strip to the whole z-plane. 
We thus get 

52 _ 



c (4) = (2?r) * r4 _ 



d 2 z 2 d 2 z 3 d 2 z 4 

-^r^r^r w*3iw)- [g,-d] 

23 



(9.11) 



21 = 1 



and by using the residual symmetry on the integration variables we are essentially reduced 
to two integrals 



C^=3(2^R^ 



[( a - tf)* Al + (a 2 + b 2 fA 2 ] 



(9.12) 



where 



At 



d 2 Z2 d 2 Z3 d 2 Z4 (\Z2 


1^3 1 


f3 2 ! 

M) 4 - 




\ z 12 1 


^34 




/3 2 

2tt 


2n 2n 2n | 


Z12\ 2 \Z34\ 2 


.1^13 1 


pl4| 


2 23 


z 24 



(9.13) 



Z 1= l 



and 



A, 



(i 2 2 2 d 2 Z 3 d 2 Z 4 ( | ^ 2 | |^3 1 1 ) 4w 



2tt 2tt 2tt 



■ X 







^24 




kl2 | 


1^14 1 


P23 


^34 



2 



^12 ^34 



(9.14) 



21=1 



These two integrals are of Dotsenko-Fateev type [|36| . By deforming the integration con- 
tours they can be transformed into products of two factors. Each factor is the sum, with 
proper trigonometric coefficients, of one-dimensional integrals of this kind 

«i 3 

/ 11 [dv^i 1 - v i) Pi ] (! - ^) 7l (l - ^ 3 ) 72 (1 - ^i^s) 73 - (9.15) 
J ° i=i 

They can be formally integrated by binomially expanding the last three factors of the 
integrand and using the fundamental integral f Q dvv a (l-v) b = T(l+a)r(l+6)/r(2+a+6). 
Therefore the two integrals Al and A2 can be reduced to the computation of (products of) 
converging series of three indices. We don't give the explicit expressions because of their 
algebraic heaviness. The method is really a straightforward generalization of the paper 



1 37]. The result can now be evaluated numerically by extrapolating the finite sums. 



The integral contributing to is evaluated with the method developed in |J8| . We 
first map the finite strip to the annulus p < \z\ < 1, where p = exp(— 277 L/R), and then we 
compute the leading contribution of the integral in the limit p — > 0. Using the symmetry 
under permutation of the four integration variables we obtain 



C< 4) = 12(2tt)^ 



lim — jy, 
p^o ln(- 



dr 4 



r-4 



dr ? 



?"3 



dr-. 



r2 a rr r 



4* Q 



2 • 



(9.16) 
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Having ordered the integration variables it is now possible to binomially expand each factor 
z Ji = r l ex P(^j7)(l — ex P(^(#z — @j) r i/ r j) 1 f° r j > I- Then we obtain a series on 12 indices 
constrained by three independent conditions from the angular integrations, each term of 
the series being a product of binomial coefficients and four radially ordered integrals of 
powers of the ?v The last integral, the one in dr 4 , gives the necessary overall volume 
divergence (|2.1J ) ln(l/p). The result is 

= 3(2tt)^R 4 -^ [(a 4 - b 4 ) 2 S 1 + (a 2 + b 2 ) 4 S 2 ] (9.17) 

where we give as an example a term contributing to S 2 



E 



u.^i... ( n i + n 2 + n 3 + | + £)(n 3 + n 5 + n 6 + \ + £)(n 2 + n 3 + n 4 + n 5 + 1 + f^) 

(9.18) 

where the symbol ^ means the following conditions on the indices ni, . . . uq, mi, . . . uiq 

ni + n 2 + n 3 = mi + m 2 + m 3 ; n\ — n 4 — n§ = mi — m 4 — m 5 

(9.19) 

n 2 + n 4 - n 6 = m 2 + m 4 - m 6 ; n 3 + n 5 + n 6 = m 3 + m 5 + m 6 . 
The coefficient g^™, is 

^-9^-g^-9>™KKKKK 

(9.20) 

with 6 n (x) = T(n + 1 — x)/ (n!T(l — x)). 

Summarizing, the fourth order correction for the running central charge is given by 

= 3(2tt)^ [(A 2 + S 2 ) k\ + S x k 2 _k\ + A x k 4 _] r 4 ~^ (9.21) 

where the numbers Ai,A 2 ,S±,S 2 depend on £- and can be computed numerically by extrap- 
olating the values of the respective finite sums. Unfortunately the above sums, especially 
the Ai ones, are slow to converge affecting therefore the precision of the extrapolated 
values. For the most favourable case § = ^ the extrapolated values are 

A 1 =49.9 ,A 2 = 1.44 

(9.22) 

Si = -20.1 ,S 2 = -1.79 . 
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For the first two coefficients A\ and A 2 we have used the VBS extrapolation method 
over the set of finite sums with 25 < N < 40, while the other coefficients £1 and S 2 
have been determined with the BST extrapolation method with convergence parameter 
u = 1 over finite sums with N < 7. The need for difference extrapolation methods is 
due to the difference in the rate of convergence of the series. With our choice the given 
extrapolated values are the most stable with respect to N. A coincise introduction to the 



above extrapolation methods can be found in |[39|| . 

The numerical integration of the TBA equations gives us the running central charge 
C(r, m/mi) with high precision and therefore the coefficients of the expansion 



C(r, m/rai) = 2 + Br 2 + 2J a 2k r k{2 ~^ 



(9.23) 



k=i 



can be determined with a standard fitting procedure. By matching the perturbative com- 
putation with the first two coefficients a 2 and 04 we can determine now the two functions 
fc_|_(m/mi), L(m/mi) as solution of the following second order algebraic system 



ci2 = 6 (2% 



0,4 




+ £1 



+ (A 2 + S 2 ) 



(9.24) 
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